Initiation of granular surface flows in a narrow channel 
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We experimentally investigate how a long granular pile confined in a narrow channel destabilizes 
when it is inclined above the angle of repose. A uniform flow then develops, which is localized at the 
free surface, ft first accelerates before reaching a steady uniform regime. During this process, an 
apparent erosion is observed and the thickness of the flowing layer increases. We precisely study the 
evolution of the vertical velocity profile in this transient regime. The measurements are compared 
with the prediction of a visco-plastic model [P. Jop, Y. Forterre and O. Pouliquen, Nature 441, 727 
(2006)]. 



A characteristic of dry granular materials is that they 
can behave like a solid or a liquid. A typical situation 
is obtained when an avalanche is triggered at the sur- 
face of a pile. In this case, grains start to move at the 
free surface, accelerate and put into motion other grains 
initially static. Understanding how the flowing part in- 
teracts with the static part has motivated many experi- 
mental works. Different configurations have been inves- 
tigated: avalanches propagating on a static layer inclined 
with respect to the horizontal jll, [3j @] i pile collapsing on 
an horizontal surface U, [ali flows in rotating drums or 
on a pile [1, S H, H, [l3jlH ■ However, the dynamics ob- 
served in these experiments is complex, since the frontier 
between flow and no-flow evolves both in space and time. 
Investigating the erosion process in a uniform avalanche, 
where the flow/no-flow interface varies only in time, is 
one of the motivation of this study. 

From a theoretical point of view, different approaches 
have been proposed to describe the solid-liquid transi- 
tion. A whole class of models is based on depth averaged 
equations and writes the mass and momentum equation 
for the flowing layer and the static pile \^ . In this frame- 
work, an additional closure equation has to be proposed 
to describe the evolution of the interface [l^ [3 Il5| . 
A second approach considers the granular material like 
a mixture of a liquid and of a solid phase and writes 
an empirical equation for the liquid phase fraction [l6l|. 
This approach captures some non trivial features of the 
transition between static and flowing regions [16]. 

Recently, it has been shown that for some conflgura- 
tions with sidewalls, the localization of the granular flows 
on top of a pile is simply related to the non uniform 
distribution of stresses. Due to the lateral friction, the 
ratio between shear stress and normal stress decreases 
when going deeper in the pile. At a critical depth, it 
reaches the yield threshold and the material stops. Us- 
ing a visco-plastic rheological model [l^jliBl' quantitative 
predictions have been obtained for steady uniform flows. 
One can then wonder if unsteady avalanches, where ero- 
sion is observed, can be captured by the same approach. 
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To study initiation of flow, we design an experimen- 
tal set-up where a long pile conflned in a narrow chan- 
nel is suddenly destabilized above the angle of repose. 
This configuration allows to create an accelerating sur- 
face flow, where the interface between flow and no-flow 
regions is flat and evolves in time only. The experimen- 
tal set up is presented in Fig. [H It consists in a long 
channel with glass sidewalls and a rough bottom to pre- 
vent the whole pile from flowing. We use glass beads 
d = 0.53 ± 0.05 mm in diameter, and density ps = 2450 
kg/m'^. In order to study a quasi-2D system we choose 
a narrow channel W = 19d (1 cm). The velocity of the 
grains measured through the glass sidewall is then rep- 
resentative of the bulk behaviour ,17]. A long L-shaped 
wedge closes the channel and defines a rectangular box 
110 cm long and 10 cm high. The wedge can rotate 
around its upper tip. At the beginning of an experiment, 
the channel is empty and the wedge is in the bottom 
position. The whole set up is inclined at the desired in- 
clination 9. The beads are then poured below the wedge 
from the top. To have reproducible initial conditions we 
gently tap several times on the set-up. Thanks to this 
set-up, a long static pile can be created at any arbitrary 
inclination [Fig. [TJa)], which allows to study initiation of 
flow above the angle of spontaneous avalanches. 

At t = 0, we release the heap by rapidly lifting the 
wedge. The time evolution of the pile is sketched in 
Fig. [TJ The lower part of the pile collapses (black ar- 
row) and creates an uphill front (white arrow). At the 
same time a uniform layer of grain starts to flow in the 
middle of the pile. Due to the downhill flow, the height of 
the rear of the pile decreases, leading to a downhill front 
(white arrow). In between these two fronts, grains accel- 
erate but the flow properties remain uniform along the 
X-axis. All the measurements presented in the following 
are done in this uniform central region (circles in Fig. [T]) 
before the two fronts arrive. The measuring time allowed 
in this uniform state varies between 1 s and 2.5 s. 

The motion of the grains in this central region is 
recorded from the side using a fast camera (500-1500 fps). 
From the movies we compute the instantaneous velocity 
proflles using a PIV method (direct correlations). The 
figure [2l[a) shows a typical time evolution of the velocity 
profile. Each profile V{z, t) is obtained by averaging in 



2 



(a) 



Rising 
system 




\ Height = 1 



= 10 cm 



Width = 1 cm 



Fast camera 
■ Rough bottom 




FIG. 1: (a) Sketch of the experimental set-up. (b) and (c) 
Sketch of the pile evolution when the wedge is suddenly re- 
leased. Measurements are carried out in the circled zone, 
before the two fronts (white arrows) reach it. The dashed 
hne represents the moving interface between the flowing layer 
and the static pile. 



the x-direction over the recorded window and averaged 
over 10 different runs carried out in the same conditions. 
Here, the origin of the z-axis corresponds to the position 
of the free surface and takes into account the slight dila- 
tancy observed during the dynamics (from 1 to 2 particle 
diameters). The first observation in Fig.[2l^a) is that the 
velocity profile evolves toward a steady shape, shown by 
the accumulation of the curves. This saturated state cor- 
responds to the steady uniform flow, which develops at 
the same angle when the beads are continuously supplied 
[TtI . The second observation is that the flowing layer get 
thicker as it accelerates, meaning that both acceleration 
and erosion take place simultaneously. Finally one ob- 
serves that the velocity profile is first exponential, but 
changes shape when the layer accelerates. Still, an expo- 
nential tail is observed deep in the pile as shown in the 
lin-log plot [inset of Fig.[2l^a)]. 

We have systematically studied the dynamics by 
performing experiments at different inclination angles 
(26.1° < 6* < 32.15°). The time evolution of the free 
surface velocity Vs{t) is shown in Fig. [S^a) for different 
angles (solid lines). The velocity increases and saturates 
when the steady regime is reached. When increasing the 
inclination, the acceleration is higher but it takes longer 
to reach the steady state. The same behaviour is ob- 
served for the flow rate Q obtained by integration of 
the velocity profile over the depth: Q{t) increases more 
rapidly for high angles but the steady regime establishes 
later [Fig.[3tb)]. 

To analyse the dynamics of the erosion process, we 
would like to define a flowing thickness. However, due 
to the exponential tails, there is no position correspond- 
ing to a strictly zero velocity. We have therefore de- 
fine two lengthscales characterizing the velocity profile: 
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FIG. 2: (a) Evolution of the velocity profile at the wall V{z) 
for Q — 28.0°. The time separating two curves is At = 30 
ms during 1.8 s. Inset: semi-logarithmic scale; (b) and (c) 
Comparison between experiments (symbols) and numerical 
simulation (dashed lines) for the velocity profiles: (b) 6 — 
26.1° at tj yj djg — 1.2, 15.1 and 166.1 (Inset: same data in 
lin-log), (c) e = 32.15° at tj y/dfg = 2.3, 7.5, 24.2, 77.8 and 
175 (Inset: same data in lin-log). 



A the characteristic length of the exponential tail and 
iji) (t) — Q/Vs the mean flow depth. With this defini- 
tion of (h), a fully exponential profile would give (h) = A, 
whereas a fully linear profile over a thickness H would 
give {h) — H/2. Comparison between (h) and A then 
gives information about how the profile differs from an 
exponential. The time evolution of (h) and A are plotted 
in Fig. m The first striking observation is that the mean 
flow thickness (h) starts from a non-zero value [Fig.UJa)]. 
This means that, as soon as the avalanche starts, a finite 
thickness instantaneously starts flowing. Then the thick- 
ness increases and eventually reaches its steady regime 
value. At the same time, the length A of the exponen- 
tial tail slightly decreases in time, but essentially remains 
constant equal to 2 particle diameters, independently of 
the inclination [Fig. [IJb)] . This value for A is consistent 
with other measurements made in steady flows 0,3 or in 
transient avalanches near the angle of repose |l9| . How- 
ever, contrary to the spontaneous avalanches studied by 
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FIG. 3: (a) Time evolution of the free surface velocity Vs for 
6 = 26.1°, 27.03°, 28.0°, 28.85°, 29.5° and 31.2°. (b) Time 
evolution of the flow rate Q; the solid lines correspond to 
measurements and the dashed lines to simulations. 



Courrech du Pont et aZ. [T^l, the velocity profile in our 
experiments conducted above the angle of avalanche is 
not always exponential. This is shown inset of Fig. |D^b) 
where the ratio (h) /A is plotted for different inclinations. 
This ratio should be equal to 1 for a pure exponential 
profile. At low inclination, the velocity profile remains 
close to an exponential during the whole process, as in 
[Tot . However, for higher inclinations, the velocity pro- 
file differs from exponential as soon as the flow starts 
(Fig. OJc) and the ratio {h) / A increases during the flow 
acceleration. 




100 150 200 



100 150 200 



FIG. 4: Time evolution of the average thickness (h) — Q/Va 
(a) and of the length of the exponential tail A (b) for 6 = 26.1° 
(★), 28.0° (A) and 31.2° (o). Experiments (sohd lines), nu- 
merical simulations (dashed lines). Inset of (b): ratio (h) /A. 



Our measurements then show that the velocity pro- 
file follows a non trivial dynamics where both accelera- 
tions and erosion are present. In the following, we try 
to compare these measurements with the prediction of 
the constitutive law proposed in [l^ . We know that the 
rheology quantitatively described the assymptotic steady 



regime [17|, and we would like to test if it can also cap- 
ture the flow initiation. The rheology assumes that the 
granular material behaves like an incompressible fluid of 
density ($ being the volume fraction equal to 0.6) 
with a visco-plastic constitutive law. The relation be- 
tween the stress tensor aij and the shear rate tensor jij 
is given by: 



with 



= |.| lij, (1) 



where P represents an isotropic pressure, I7I represents 
the second invariant of the shear rate tensor j-yj — 

and p{I) is the friction coefficien 0, [l^, HH, 

which depends on the incrtial number / = as fol- 

lows: 



KI) = Ms + 



p.2 - ^is 



(2) 



In the configuration studied here, the velocity is par- 
allel to the free surface and depends on z and y, v — 
V{y, z, t)e^, and the pressure is given by P = Ps^gz cosO. 
To simulate our experiments, we then have to solve the 
Cauchy equation: 

Ps^dV/dt = ps^gsinO + da^z/dz + duxy/dy, 

with the following boundary conditions: free-stress at 
the free surface, no-slip condition at the bottom and a 
Coulomb friction on sidewalls. To quantitatively com- 
pare the prediction with the experiments, we choose 
p, = tan(20.9°), p2 = tan(32.76°) and Jq = 0.279 as 
used in pL7|. We also need the value of the wall friction 
coefficient pw We choose pw = tan(13.1°), which gives 
the best agreement for the mean flow rate in the steady 
uniform regimes based on data given in (l7| . With this 
choice, no free parameter remains to quantitatively pre- 
dict the dynamics of our avalanches. 

The simulation is carried out by imposing a zero initial 
velocity and by letting the system evolve under gravity 
(Fig. [5]) . At t = a finite layer of material starts flow- 
ing. This layer then accelerates and the flow thickness 
increases in qualitative agreement with the experiments. 
In order to quantitatively compare the simulation with 
the experiments, we extract from the simulation the ve- 
locity profile at the sidewalls and compare with the ex- 
perimental data. The prediction are plotted in Figs.[2{b) 
and die) (dotted lines) for 9 = 26.1° and 9 = 32.15°. We 
see that the agreement for the evolution of the velocity 
profile is good at high inclination, but becomes less ac- 
curate at low inclination. In the model the velocity goes 
to zero at a finite depth, whereas in the experiments the 
profiles are followed by an exponential tail. At low angles 
the experimental profiles are mainly exponential and the 
model does not capture the correct shape. 

In order to more precisely check if the dynamics is cor- 
rectly predicted by the model, we extract from the simu- 
lations the same quantities computed in the experiments. 
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FIG. 5: Velocity profiles from the simulation for 9 — 26.1° 
at different times t/y^d/g = 1.5, 13.4 and 200 (stationary 
regime). The solid lines represent the frontier between the 
flowing layer and the static bottom for the three profiles. 



namely the free surface velocity at the wall Vg (t) , the flow 
rate at the wall Q{t) and the mean thickness (h) (t). The 
results are the dotted lines in Fig.[3]and Fig.lH We find a 
good agreement for the free surface velocity and the flow 
rate, the time evolution being quantitatively captured by 
the model for all inclinations. The agreement is less ac- 
curate for the mean thickness (h) (t). As shown in Fig.lH 



the model predicts that (h) (t) is initially non zero and 
increases with time in a similar way the experiments do. 
However, quantitatively the prediction is not very accu- 
rate. The major reason for this discrepancy is that the 
parameter (h) (t) is sensitive to the exponential creeping 
tail, which is not captured by the model. 

In conclusion, we have shown that in a narrow channel 
configuration, for which wall effects are predominant, the 
observed erosion process is captured by a simple visco- 
plastic model. The inertia and the existence of a thresh- 
old are sufficient to create an "erosion- like" phenomenon. 
The frontier between flow and no flow regions is simply 
driven by the time evolution of the internal stresses. The 
model quantitatively predicts the time dynamics show- 
ing the relevance of the proposed rheology for unsteady 
flows. On the other hand, serious limits exist. First, 
the model fails in precisely predicting the velocity proflle 
shapes, mainly because the creeping exponential tail is 
not captured. The discrepancy is even more important 
at low inclinations, close to the angle of repose. The sec- 
ond limit concerns the lack of hysteresis in the model. 
Our study is restricted to avalanches triggered above the 
angle of repose. To capture natural avalanches observed 
for example when slowly inclining a pile, or to describe 
intermittent flow regime, hysteresis has to be taken into 
account in the model. These limits shows that further 
developments are needed for a more complex rheology, 
which could capture both the quasi-static limit and the 
hysteretic nature of granular flow threshold. 
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